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Abstract 

Melonic graphs constitute the family of graphs arising at leading order in the 1/N expansion of 
tensor models. They were shown to lead to a continuum phase, reminiscent of branched polymers. We 
show here that they are in fact precisely branched polymers, that is, they possess Hausdorff dimension 
2 and spectral dimension 4/3. 

1 Introduction 

Tensor models [1] have recently emerged as a promising conceptual and computational playground for the 
modern theoretical physicist. Drawing ingredients and impetuses from several topics in pure and applied 
mathematics, they provide a theory of random higher-dimensional spaces. More precisely, tensor models 
are theories of tensors in the fashion of matrix models [2], of which they may be viewed as a superset. 

In particular tensor models have been used in non-perturbative quantum gravity [3, 4, 5, 6, 7, 8, 9, 10, 11, 
12, 13]. In this context, attention concentrated mostly on building models well-adapted to the problem 
of interest, namely, a quantum theory of gravity. However, within these models, computations have been 
notoriously difficult to perform. There are two reasons. Firstly, most models initially studied were highly 
complex, attempting to incorporate all the features developed over many years of research in the area of 
discrete quantum gravity [14, 15, 16, 17]. To date, this has precluded any analysis beyond the 1st order in 
perturbation theory [18, 19, 20]. Secondly, for a long time, no generic tool existed in practice, with which 
to extract information about their generic features at large orders in perturbation theory. 

This situation changed with the advent of what one might call modern tensor models [21]. Initially striking 
for their apparent simplicity, on closer inspection, they possess a rich and nuanced structure. In particular, 
they allow for the development of a 1/iV-expansion [22, 23, 24, 25, 26], that is, they come equipped with a 
parameter N that sensibly organizes the Feynman graphs according to its inverse powers. In the large-TV 
limit, one finds that only a certain subclass of graphs survive and contribute to the free energy, the melonic 
graphs. It is on this melonic sector that most research has concentrated, resulting in the extraction of 
(multi-) critical exponents [27]. 

These results have excited a period of high productivity, giving credence to the claim that the framework 
facilitates explicit calculations involving large numbers of highly refined graphs. The iro-class of models 
has been extended to involve matter degrees of freedom of both Ising/Potts [28] and hard dimer types 
[29, 30], and an extensive analysis of their critical behaviors in the melonic sector has been performed. 
A class of dually weighted models [31] has been identified with the same melonic analysis conducted. 
The quantum symmetries have been formally identified [32, 33, 34, 35] and catalogued at all orders and 
the resulting generators conform to a higher-dimensional Virasoro-esque Lie algebra. Indeed, there are 
even universality results for the critical behavior [36]. This foundational work has influenced another 
concurrent line of research on tensor field theories, that is, tensor models with modified propagators that 
induce a renormalization group flow [37, 38, 39]. These tools may be expected to be of use in analyzing 
the gravity-inspired models mentioned earlier [40, 41]. 

Intriguingly, the critical exponents extracted for the melonic sector in the plain, matter-coupled and 
dually weighted models are identical to those for branched polymers [42, 43], which are a class of trees. 
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Furthermore, both melonic graphs and branched polymers have spherical topology. However, while mclonic 
graphs are reminiscent of branched polymers, these two structures are not naively identical. As graphs, 
they have markedly different connectivity and thus, need not exhibit the same intrinsic physical properties. 

Fortunately, branched polymers are very well-studied. After all, they have been seen to characterize phases 
in both euclidean and causal dynamical triangulations [44, 45, 46]. In particular, there is solid information 
on their Hausdorff [42] and spectral [47] dimensions, which are dn = 2 and ds = 4/3, respectively. In 
this paper, we show that melonic graphs have the very same Hausdorff and spectral dimensions. What 
is more, we claim that this conclusively identifies the class of melonic graphs as branched polymers. One 
might wonder why these are sufficient criteria. This stems from the fact that both the Hausdorff and 
spectral dimensions have a deep physical significance. The Hausdorff dimension governs the scaling of 
the volume with respect to the geodesic distance, while the spectral dimension is the effective dimension 
experienced by a diffusion process. Thus, both observables are intrinsically linked to physical processes on 
the structures in question. 

Fortunately, we may co-opt the large body of literature on branched polymers for our own purposes. In 
other words, arguments translate from that case into this more general one. Having said that, certain 
details change upon generalization and we make efforts to highlight these clearly as we progress through 
the reasoning and calculations. 

Two papers, in particular, have greatly influenced this work. The first [48] deals with the Hausdorff 
dimension of (2-dimensional) stack triangulations. These are particular triangulations of the 2-sphere 
that are in correspondence with branched polymers. The authors show in particular that the Hausdorff 
dimension of stack triangulations equals 2. Thus, they faced an analogous problem to ours; to show that a 
class of objects that were, at the outset, merely in correspondence with branched polymers, truly possessed 
the same physical property. As a result, it was an invaluable guide through this murky terrain. The second 
[47] deals with the spectral dimension of branched polymers directly. Although, the onus of generalization 
fell directly on us, this work provides both the solid background and numerous helpful insights through 
various technical challenges. 

This paper is organized as follows. We introduce in Section 2 (with details and references in Appendix 
A) the class of melonic graphs within the tensor model setting, before proceeding to set up various cor- 
respondences to other useful classes of objects, in particular, to a class of simplicial D-balls and a class 
of (D + 1 )-ary trees. In Section 3, we construct metrics on these two classes of objects, concentrating on 
vertex depth, that is, the distance from a generic vertex to a distinguished root vertex. We begin Section 
4 by showing that the depth of a randomly chosen vertex in the D-ball is, in the infinitely refined limit, 
a fixed rescaling of the depth of the corresponding vertex in the associated (D + I)-ary. This result is 
important for the subsequent argument (augmented by further explanation in Appendix C) pertaining to 
the Hausdorff dimension of the melonic Z?-balls. Finally, Section 5 deals with the spectral dimension of 
melonic graphs. 

2 Tensor models and melonic graphs 

As one might imagine, tensor models are built from tensors: collections of N D complex numbers, where N 
is the size of the tensor and D is the dimension, that is, the number of indices. The models themselves are a 
specific class of perturbed-Gaussian probability measures that are independent and identically distributed 
(iid) across the components of the tensor. In order to extract information about these measures, one 
should evaluate their partition functions and moments. Let us say a word about the partition function. 
For small values of the perturbation parameter, one may expand the measure as a Taylor series 1 . After 
evaluating the Gaussian integrals using Wick contraction, one has a sum of terms labeled by Feynman 
graphs. The Feynman graphs of the iid measure are bipartite edge colored graphs, see Appendix A for 
more details. All the vertices of the graphs have valence D + 1 and all the edges have a color 0, I up 
to D such that the D + 1 edges incident at a vertex have distinct colors. It emerges that these terms 
may be organised according to the their power of l/N (which is alway non-negative), hence the name 
l/W-expansion. The iid measure encodes a uniform distribution over the graphs at any fixed power of 
l/N. Moreover, as N f oo, only one subset of graphs survives; those for which the power of l/N is zero. 
These graphs are known as melonic graphs, a name which stems from their distinctive shape. In fact, 
we shall analyse the properties of a related set of graphs known as rooted melonic graphs. These are 
melonic graph with one edge of color cut and they yield the 2-point function at leading order in the 

lr The ensuing perturbation series is not summable, but it sometimes is Borel summable [36] 
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1 /TV-expansion. 



In short, we describe the structure of rooted melonic graphs below, along with their relation to rooted 
(D + l)-ary trees, melonic D-balls and stack D-spheres. We refer the reader to Appendix A and references 
cited therein for details. 

2.1 Connected 2-point function and rooted melonic graphs 

Rooted melonic graphs have slightly more structure than their non-rooted counterparts, are more versatile 
and hence, prove easier on the whole to work with. They are defined in an iterative manner. 

The fundamental building blocks of any rooted melonic graph are the elementary melons. Such a melon 
consists of two vertices connected by D edges. Both vertices have one external edge. Obviously, both 
external edges possess the same color, say i (and one refers to such an object as an elementary melon of 
color i). An elementary melon has two distinguished features: i) an external edge of color i incident to 
the white vertex, which is known as the root edge; ii) D + 1 edges incident at the black vertex, which 
are known as active edges, having distinct colors from {0, 1, . . . D}. The iterative definition proceeds as 
follows: 

p = 1: There is a unique rooted melonic graph with two vertices. It is illustrated in the bottom left of 
Figure 1 and is the elementary melon of color 0. 




Figure 1: An elementary melon of color 2 inserted along the active edge of color 2 (for D = 3). The 
active edges are drawn using full lines. 

p = 2: There are D + 1 melonic graphs with four vertices. One obtains them from the graph at p = 1 
by replacing an active edge of a given color by an elementary melon of the same color (as shown in 
Figure 1). 

p = k: One obtains these graphs from those at p = k — 1 by replacing some active edge by an elementary 
melon of the appropriate color. 

2.2 The many guises of rooted melonic graphs 

As mentioned earlier, the abstract structure of rooted melonic graphs coincides with that of several other 
objects, which we shall describe presently. 

2.2.1 Colored rooted (D+l)— ary trees 

There is a simple bijection between the set of rooted melonic graphs and colored rooted (D + l)-ary 
trees. The fundamental building blocks of any colored rooted (D + l)-ary tree are the elementary 
vertices. An elementary vertex of color i is (D + 2)-valent with two distinguished features: i) a root 
edge of color i; ii) D + l active leaves each with a distinct color from {0, . . . , D}. These correspond to 
the root edge and the active edges of the elementary melon, respectively. Since this class of trees is also 
constructed in an iterative manner, the map is self-evident: 

p = 1: There is a unique colored rooted (D + l)-ary tree with a single elementary vertex. This is the 
elementary vertex of color 0. It is illustrated in the bottom left of Figure 2. 
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Figure 2: The elementary vertex of color 2 replacing a leaf of color 2 (for D = 3). 

p = 2: There are D + 1 such trees with two elementary vertices. One obtains them from the tree at p = 1 
by replacing a leaf of a given color with an elementary vertex of the same color (as shown in 
Figure 2). 

p = k: One obtains these trees from those at p = k — 1 by replacing a leaf with an elementary vertex of 
the same color. 



2.2.2 Colored simplicial D balls 

The description of how any given (D + l)-colored graph is dual, in a precise topological sense, to a unique 
Z?-dimensional abstract simplicial pseudomanifold is provided in all its detail in many contexts [49, 50] 
(see [51, 52] for a modern description). We give a heuristic summary here. 

Consider first a closed (D + l)-colored graph. One identifies the set gw— **) of all maximally connected 
subgraphs with k distinct colors . . . , ik} drawn from {0, . . . , D}. These are known as the k bubbles 
of species (ii, . . . , ik). Note that these bubbles have a nested structure in that /c-bubbles lie nested within 
(fc + l)-bubbles, which in turn lie nested within (fc + 2)-bubbles and so on. 

The dual map associates a (D — fc)-simplex to each fc-bubble and the nested structure of the bubbles 
encodes the gluing relations of the simplices to form a unique simplicial complex. These (D — /c)-simplices 
inherit the coloring of their corresponding fc-bubble. Note that the vertices (0-simplices) of the dual 
simplicial complex are one to one to the subgraphs with D colors. The vertices are thus colored by D 
colors {ii, . . . id}- Alternatively, one can color them by an unique color: the complementary color of their 
dual subgraph, {0, . . . D} \ {ii, . . .in}- 

Consider cutting a closed (D + l)-colored graph along one edge. This results in an open graph whose 
dual is a simplicial complex with boundary. This boundary is a {D — l)-sphere constructed from two 
(D — l)-simplices. 

Melonic graphs are dual to simplicial D-spheres. Rooted melonic graphs, which are melonic graphs with 
one edge cut, are dual to simplicial D-balls with the boundary mentioned above. For the want of a 
better name, we shall call them melonic D balls. One can define them iteratively. The fundamental 
building blocks are the elementary melonic D balls. These consist of two £)-simplices sharing D of 
their (D — l)-simplices. There are two more (D — l)-simplices forming the boundary [D — l)-sphere. In 
the manner outlined above, the two D-simplices are dual to the two vertices of the elementary melon, 
while the (D — l)-simplices are dual to the edges. Thus, the (D — l)-simplices inherit a single color. An 
elementary melonic D-ball of color i has two distinguished features: i) an external root (D — l)-simplex 
of color i; ii) D + 1 active (D — 1)— simplices (one of which is on the boundary), each with a distinct 
color. The iterative definition proceeds as follows: 

p = 1: There is a unique D-ball comprised of two D-simplices. It is the elementary melonic Z?-ball of 
color 0. It is illustrated in the bottom left of Figure 3. 
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Figure 3: The melonic D-ball at p = 2 (for D = 3), obtained by adding an elementary D-ball of color 2. 

p = 2: There are D + 1 melonic D-balls with four D-simplices. One obtains them from the melonic D- 
ball at p = 1 by adding an elementary melonic D-ball of a given color (shown in Figure 3). More 
precisely, one splits the melonic D-ball arising at p = 1 along an active (D — l)-simplex (or one 
selects the active boundary (D — l)-simplex). One then glues an elementary melonic D-ball of the 
appropriate color along the split (or simply on the boundary (D — l)-simplex). 

p = k: One obtains them from those at p = k — 1 by adding an elementary melonic D-ball at some active 
(D — l)-simplex. 

2.2.3 Colored stack simplicial D-spheres 

There is also a slightly more convoluted map between the set of colored rooted (D + l)-ary trees and col- 
ored stack simplicial D— spheres (abbreviated here to stack spheres). Stack spheres are triangulations 
of the D-sphere. The fundamental building blocks are the elementary stack spheres. Such an object 
is comprised of the D + 2 D-simplices forming the boundary of a (D + l)-simplex. An elementary stack 
sphere of color i has two distinguishing features: i) it has a root D simplex of color i; ii) D + l distinctly 
colored active D— simplices drawn from {0, . . . , D}. Moreover, there is a single vertex which is shared 
by all D + 1 active D-simplices, which we shall refer to as its active vertex. An example is illustrated in 
Figure 4. 

From this coloring of the £)-simplices, one can color the vertices in the following manner. The root D— 
simplex of color i contains D + 1 of the D + 2 vertices. Consider such a vertex. There is a unique active 
D-simplex, within which it is not contained. It is labeled by the color of that D-simplcx. The final vertex 
(the active vertex) is not labeled by a unique color but rather by all colors. Better said, the color of this 
vertex depends on the active D-simplex, within which one is considers it, and one labels it by the color of 
this D-simplex. One notes that for any given D-simplex in the elementary stack sphere, its D + 1 vertices 
are distinctly colored. This coloring procedure is illustrated in Figure 5. 




Figure 5: Coloring the vertices of the elemen- 
Figurc 4: The unique stack D-sphere at p = 1 tary stack sphere (for D = 2). Only the active 
(for D = 2). 2-simplices are drawn. 

p = 1: There is a the unique stack sphere comprising of D + 2 D-simpliccs. It is the elementary stack 
sphere of color 0. 
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p = 2: There are D + 1 stack spheres with 2D + 2 D-simplices. One obtains them from the stack sphere 
at p = 1 by constructing the connected sum of this stack sphere with an elementary stack sphere of 
some color. To be more precise, one takes the stack sphere at p = 1 and excises an active _D-simplex 
of some color, say i. Having done that, one removes the root £)-simplex from an elementary stack 
sphere of color i. One glues the resulting D-balls together by identifying their boundaries, such that 
the colors of their respective boundary vertices match and the result is a D-sphere. Note that this 
is akin to performing a 1 — > {D + 1) Pachner move on the active D-simplex with precisely inherited 
color information. This process is drawn in Figure 6. 



Pachner move 



Figure 6: The 1 —>(£> + 1) Pachner move performed on the active D-simplex of color 2 (for D = 2). 

p = k: One obtains these stack triangulations (which have Dk + 2 D-simplices) from those at p = k — 1 
by performing a 1 — > (D + 1) Pachner move on an active -D-simplex. 

2.2.4 Words on trees 

In principle, the vertices of a colored rooted (D + 1)- 
ary tree are of two types. There are univalent ver- 
tices, at which the leaves and root edge of the tree 
are incident. These vertices are disregarded from now 
on. Rather, one concentrates on the (D + 2)-valent 
vertices, to which we shall continue to refer as ele- 
mentary vertices. Each tree also has a distinguished 
elementary vertex known as its root vertex. 
Every elementary vertex of a colored rooted (D + 1) 
ary tree has an associated word composed from the 
alphabet Hd = {0, 1, . . . , D}. For a given vertex, its 
word is constructed by listing (left to right) all the col- 
ors one encounters on the branches, when going from 
the root vertex to the vertex in question. Thus, the 
word associated to the root vertex is ( ). An exam- 
ple of a generic word is given by (10132120312). This 
word is the canonical label discussed in [27]. It will 
serve our purposes better to attach to the vertices, 
this word supplemented by an initial letter 0, which 
signifies the initial root edge of color 0. Thus, the la- 
bel of the root vertex becomes (0; ), while the sample 
word given a moment ago becomes (0; 10132120312). 
This example is drawn in Figure 7. 



Figure 7: The words associated to some 
vertices of a colored rooted (D + l)-ary 
tree. 

2.2.5 More remarks on vertices in trees, balls and stack spheres 

At order p = k, a rooted melonic graph has 2k internal vertices, the associated colored rooted (D + l)-ary 
tree has k elementary vertices, the associated melonic Z?-ball has k internal vertices, while the associated 
stack D-sphere has k active vertices. We shall discard any adjectives from now on and refer to these nodes 
simply as the vertices. 






(0; 10132120312) 

(0; 1013212031) 

(0; 101321203) 
(0: 10132120) 
(0; 1013212) 

(0; 101321) 

(0; 10132) 

(0; 1013) 
(0; 101) 

(0; 10) 
(Oil) 

(0; ) 
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There is a clearcut correspondence between the vertices of the melonic D-ball and the vertices of the 
(D + l)-ary tree. Both the D-ball and the (D + l)-ary tree are generated iteratively with a single vertex 
added at each iteration. Moreover, both have a distinguished initial vertex called the root vertex. As 
mentioned already, in the (D + l)-ary tree, it is the vertex labeled by the word (0; ). In the U-ball, it is 
the corresponding vertex, which is of color and denoted by vq ■ 

In the same fashion, there is a correspondence between the vertices of the (D + l)-ary tree and the vertices 
of the associated stack Z?-sphere. For completeness, there is a two-one correspondence between the vertices 
of the (D + l)-ary tree and the associated rooted melonic graph; each vertex of the tree is associated to 
an elementary melon, which in turn has two vertices. 

3 Distance and depth 

3.1 Distance 

For any connected graph, there is an elementary definition for the graph distance d(v s , v t ) between two 
vertices v s and Vt as the minimal number edges in any contiguous path journeying from v s to Vt- 

In this case, there are four related graphical structures: a rooted melonic graph, a colored rooted (D + 1) 
ary tree, a melonic D-ball and a stack Z?-sphere. Although their vertices correspond in the manner 
detailed above, they have different graph connectivities. Therefore, the distance between a pair of vertices 
in the (D + 1)— ary tree differs in general from the distance between their corresponding vertices in the 
melonic £)-ball and so forth. 

3.2 Depth 

In Section 2.2.4, we distinguished the vertex at the base of a (D + l)-ary tree as the root vertex. There are 
corresponding root vertices in both the associated melonic D-ball and stack D-sphere. We are interested 
in calculating the graph distance of a vertex from the root vertex. 

The term tree depth dx shall refer to the graph distance of a vertex from the root in the (D + l)-ary 
tree. In passing, this is rather simple to calculate since there is a unique path from the root to any 
given vertex. As mentioned above, to each vertex of a tree, there is an associated word, for example, 
(0; 10132120312). The letters to the right of the semi-colon are exactly the edges of the tree joining the 
vertex to the root. Thus, one may calculate the tree depth of this vertex by counting the number of letters 
after the semi-colon; d T (0; 10132120312) = 11. 

The term depth d shall refer to the graph distance of a vertex from the root in the melonic Z?-ball, while 
the the term stack depth shall refer to the graph distance of a vertex from the root in the stack D-sphcrc. 

Having already demonstrated how to calculate the tree depth, we shall detail here the corresponding 
procedure for the depth in the melonic Z?-ball. For the stack depth, we refer the reader to [48]. There 
are four parts. First, we give just a blind statement of the mechanism by which one constructs the depth. 
Second, we provide a worked example. Third, we explain the mechanism in more detail. Fourth, we 
describe a shorter method to calculating the depth for any given word. 

3.2.1 Construction 

As stated many times already, to each vertex of the melonic D-ball, there is a corresponding vertex in the 
associated (D + l)-ary tree. Thus, each vertex in the melonic D-ball is labeled by a word, for example, 
(0; 10132120312). It emerges that one may calculate the depth of a vertex in the melonic D-ball using the 
associated word. The key object is the distance array associated to that vertex, an array with D + 1 
entries that may be constructed directly from its word. It is from this distance array that the depth of the 
vertex in question is extracted. 

To be clear at the outset, the aim of the game is to construct the depth of a vertex in a melonic D-ball. 
However, one first considers the associated (D + l)-ary tree and the corresponding the vertex v n . Assume 
that for this vertex, the tree depth is: dr(v n ) = n. Thus, the associated word to this vertex is of the form: 
w n = (0; U\ . . . u n ), where itj are letters drawn from the alphabet S/j+i = {0, • ■ • , D}. In the tree, there 
is a unique path, of length n, from the root vertex to v n . Let us denote the root vertex by vq. This path 
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comprises of a sequence of n + 1 vertices marked Vi, such that < i < n. The word labelling the root 
vertex w is w o — (0; ), while the word labelling Vi is just w n with the sub- word u i+ i . . . u n removed, that 
is to, = (0; U\ . . . Uj). Moreover, given that one has a sequence of vertices in the (D + 1)— ary tree, one also 
has a corresponding sequence of vertices in the melonic D—bell which inherit the names Uj. 

One constructs the distance array associated to Vi iteratively from the words Wj, with < j < i. At the 
outset, one is given that the distance array associated to v is Oa(uo) = (0, 1, . . . , 1), that is, zero followed 
by D ones. Now, let us assume that the distance array associated to is 0o(«j) = (mo, . . . , m D ) and that 
Ui+i = j. Then, the distance array associated to v i+ i is 

da(v i+1 ) = (m , mj_i, min(m fc ) + 1 , m j+1 , m D ) . (1) 

In this manner, one can construct the distance array associated to all the Vi in the sequence. 

Say that that = j and that for Mj, the distance array is 5a(uj) = (mo, . . . ,m,j, . . . ,m D ). The depth of 
Vi in the melonic D-ball is d(vi) — mj. 

3.2.2 Worked example 

We consider our favorite example, a vertex v labelled by the word w — (0; 10132120312). 
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Table 1: Calculation of distance array and depth. 



3.2.3 Explanation 

Consider the sequence of vertices Vi (0 < i < n) in the melonic D-ball such that the word associated 
to Vi is Wi = (0; iti . . . Ui). This sequence is generated in the following manner. One starts with initial 
elementary D-ball of color and adds in sequence elementary D-balls of colors m, for < i < n. In 
terms of the 1-skeleton of the melonic .D-ball, adding the elementary D-ball of color Ui entails adding 
the vertex m, (of color u 4 ) and joining it to D vertices (already present) labelled by distinct colors from 
the set {0, 1, . . . D} \ {ui}. This action distinguishes D + 1 vertices in the melonic D-ball. One defines 
the distance array associated to Vi as the array with D + 1 entries comprised of the depths of these D + 1 
vertices, ordered according to their color. 

Consider the initial elementary D-ball of color 0. This has a unique internal vertex of color 0, which is 
the root vertex v . It also has D boundary vertices labelled by distinct colors from the set {1, 2, . . . D}. 
The word associated to the root vertex is wq — (0; ), while we stated that the distance array is Oa(t>o) = 
(0, 1, . . . , 1). This encodes that the fact that the root vertex (of color 0) is of depth 0, while the D boundary 
vertices are all of depth 1. Note that one knows the depth of all vertices in this melonic D-ball. 

Now, consider the point where one has added the elementary D-balls in the sequence up to some point 
i. Assume that the distance array of Vi is 5a(i>j) = (mo, . . . , mo)- According to the above definition, this 
encodes the depths of D + 1 vertices in the melonic D-ball: the vertex Vi along with the D vertices to which 
it is connected. Within this set of vertices, one denotes momentarily the vertex with depth m k by V k . 
Now, one adds the elementary D-ball of color j = Mj+i. At the level of the 1-skeleton, this corresponds 
to adding the vertex v i+1 of color j and connecting it to D vertices labelled by distinct colors from the 
set {0, 1, . . . D} \ {j}. By direct inspection, the vertex Mj+i is joined to D of the vertices 14, that is, all 
except Vj. This means that all entries except the jth entry of J>a(t>j+i) coincide with those entries in 5a(i>i). 
Moreover, any path joining v i+ i to the root vertex must pass through at least one of these vertices Vfe. 
Thus, the jth entry is (1 + the minimum of the depths of the vertices 14 with k ^ j). 
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3.2.4 Calculation via sub-words 

Despite this iterative calculation via a distance array, there is a yet more succinct method, which shall 
play an active role later. 

Let us denote by Wd+i the set of words containing all the letters of the alphabet Sd+i = {0, 1, • • • D}. 
Consider a vertex v labelled by the word w = (0; u\U2 . . . u n ). The depth of v corresponds to division of w 
into disjoint adjacent sub-words r r , comprised of letters of depth r. Thus, To = 0. Then, t% = u\ . . .u ai , 
with Mi, U2, ■ ■ ■ u ai ^ and u 01 _|_i = 0. Furthermore, r r , for r > 1, may be one of two forms: i) r r — 
V-i+i ■ ■ ■ u a r sucn that r r Wd+i but r r u ar+ i e W D+1 ; ii) r r = u (lr _ 1+1 ...u n ii V-i+i • • • u n ^ Wd+i- 
This second possibility accounts for the fact that the last subword might be incomplete. 

The depth of a vertex v with with word w = tqT\ . . . rfc, is: 

A(w) = d(v) = k . (2) 

As an example, take the (3 + l)-ary tree illustrated in Figure 7 and Table 1. We have listed the depth 
and array labels of all the vertices in the unique path from the root vertex to the vertex with label 
(0; 10132120312). The vertical dividers highlight those letters of the word at which the depth increases. 

To track the way in which the depth is updated from left to right, one notes that as long as one does 
not encounter a second letter 0, the array label remains (0, 1, ... , 1). Hence, the depth is 1. The depth 
increases to 2 at the second occurrence of the letter (that is the first occurrence of after the semi- 
column) and the array label of this vertex is (2, 1, ... 1). Then, starting from this 0, the first occurrence 
of a letter j shifts rrij in the array label from 1 to 2. A subsequent occurrence of the letter j has no effect 
on the array. At this point, as long as there are at least two entries in the array label that equal 1, the 
depth of the corresponding vertices remains 2. The depth increases to 3 only when the array has: i) a 
unique entry rrij — 1, ii) rrii — 2 for all i ^ j and Hi) the letter one encounters is j. In this instance, the 
array label becomes rrij = 3 along with rrn = 2 for i ^ j. This occurs when, starting from the second 0, 
one has encountered (at least once) all the letters of the alphabet. Note that, very importantly, the depth 
changes exactly when one encounters for the first time all the letters in the alphabet, that is the depth of 
(0; 10132) is 3. Using the same argument, the depth changes to 4 when, starting from the first letter of 
depth 3 (which in this case is 2), one has encountered again (at least once) all the letters of the alphabet, 
hence at (0; 101321203). 

The division in sub-words goes as follows: 

w = (0; 10132120312) = (0)(1)(013)(2120)(312) (3) 

Remark on comparison to stack depth: Note that the depth defined here is different from the stack 
depth (very similarly) defined in [48] for D = 2. The following is a convincing example. (In order to 
compare the two distances more readily, we consider a (2 + l)-colored model, such that the words are 
comprised of letters from the alphabet {1, 2, 3} and the first letter 1. In other words, u = (1; u% . . . u n ).) 
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In principle, a lot of the same structures arise in both definitions. Both depths are constructed from the 
words attached to the vertices of colored rooted (2 + l)-ary trees. While the depth defined here measures 
the distance of some vertex from the root vertex in the associated melonic 2-ball, the Albenque-Marckert 
depth [48] measures the distance of the vertex from the root vertex in the associated stack 2-sphere. Since 
the connectivity of a melonic 2-ball differs from that of its corresponding stack 2-sphere, one should not be 
surprised that the depths differ. Starting from the word attached to a vertex of the tree, the procedure to 
calculate the stack depth of the associated vertex in the stack 2-sphere is almost identical to that utilized 
here to calculate the depth of the associated vertex in the melonic 2-ball. As one might expect, one defines 
and updates a(n analogous) distance array and ultimately reads off the stack depth. The only difference 
is that the stack depth gets updated after one has encountered for the first time all the letters in the 
alphabet. Thus the stack depth lags behind the depth in the melonic D-ball. 
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4 Hausdorff dimension 



For a metric space X, the Hausdorff dimension dn captures how the volume of a ball scales with respect 
to its geodesic distance. More formally defined as: 

d H = inf{d > : U d (X) = 0} , (4) 

where l-Ld{X) is the d-dimensional Hausdorff measure on X, that is: 

H d (X) = inf j<5 = ^ rf : the indexed collection of balls of radius n cover X} . (5) 

i 

Obviously, for the simple example of flat D-dimensional Euclidean space: Vjj ~ r D . Thus: 

d H = D. (6) 

It is clear however, that some work must be done to extract this dimension for the class of melonic D-balls. 
One may follow the argument in [48], which successfully navigates this ground for the 2-dimensional case, 
that is, melonic 2-balls (or rather stack 2-spheres). A complete statement of their argument would be a 
rather laborious task and would, for the most part, amount to a literal restatement. There are, however, 
some points where our D-dimensional argument differs from their 2-dimensional one, summarized by the 
fact that we analyze melonic D-balls rather than stack D-spheres and thus utilize the depth rather than the 
stack depth. This motivates Lemmas 1 and 2 below. Apart from dwelling on this important technicality, 
we content ourselves with a descriptive account of the main theorem, which offers as a by-product the 
Hausdorff dimension. Some more details are given in Appendix C 

The investigation of the Hausdorff dimension for melonic D-balls requires a technical preamble. One 
denotes by 5 1 ™ a certain sum of multinomial coefficients: 

n 1 + ...n q =n 

Sn q = E | ' (7) 

ni,...n g >l H 

and sets Sq = 0. Remark at this point that S™ = 1 and — 2" — 2. In general, one has: 
Lemma 1. 

E (- 1 ) 9_r f?V"- (8) 



0<r<q 



Proof: See Appendix B. □ 

This technicality is of immediate use in something more practical below. To set up the following lemma, one 
should recall some points. Consider a vertex v in a melonic D-ball along with its word w = (0; u\ . . . u n ). 
One knows that the tree depth of v in the associated (D + l)-ary tree is simply dx(v) = n, while if 
w = t ti . . .Tfc, then the depth of v in the melonic D-ball is d(v) = A(w) = k. Generically, there is not 
a simple formula relating dr(y) and d(v), that is, the tree depth with the depth. However, one can ask 
whether their average ratio = A(w)/n approaches some value asn-> oo. Said more precisely: 

Lemma 2. Let U\, . . . ,U n be a sequence of random variables uniformly drawn from T,d+1> and denote 
w = Oui . . . u n . One has: 

l -K{w) A A , A A ! = (D + 1) J2 i-V D - r H (D + [_ r) 2 ■ ( 9 ) 

Proof: A^ 1 is the average length of the Tj for infinitely long w, that is, A^ 1 = (\ri\). Recalling the 
definition of the Tj, one notes that |tq | = 1 and |ri| is special: |ri| = (-1 + length of the sequence of letters 
that ends the first time the letter appears). The probability P n that this sequence ends after n letters is: 

Pi = P 2 = (l-P 1 )Pi, P 3 = (1-P 1 -P 2 )P 1 
Pn = {I -Pi Pn-l)Pl = (1 - -Pi - (1 - Pl)Pl - ■ 

= (l-Fr)™-^. 



(1 - Plfpl 



{1-P 1 ) n - 2 P 1 Pj 



(10) 
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Hence: 



= £( n - 1)(1 - P 1 ) n - 1 P 1 = P X E n{l - PiY = Px(l - p M- d Pi) 1 _n_p ) 

n>l n>0 ^ 1 ' 



n>l 
I -Pi 
Pi 



(11) 



= D < oo 



The finiteness of (|ri|) implies that A A 1 is the average length of Tj, for i > 2. 

One denotes by P™ the probability that after n draws one has obtained a sequence of letters with exactly 
q distinct colors. It follows that: 



<N) = 5> - i)pg- 1 Pi = ^ E ^ 



(12) 



n>l 



ri>0 



One needs a more explicit form of P q . To this end, one notes that the probability of any fixed sequence 
of results on n draws is rpjlT" ■ The number of configurations of n letters having n\ times the color i\, n-i 
times the color i 2 up to n q times the color i q is the multinomial coefficient: 



nil . . . n q \ 



As a result: 



pn _ 

q (D + l)™ 



D + l 
Q 



ni + ...rtg-n 



E 



i ni l...n q \ (D + l)" V q 



D + l 



E (-D 

0<r<ij 



r 



(13) 



(14) 



with the help of Lemma 1. There are two checks one needs to run on this formula. First, one checks that 
the probabilities are normalized: 



D+l 

En" 

q=0 



-, -D+l 

i V 

+ 1)" ^ 



D + l 
q 



q=0 

D+l r D + l 



(D + l) 

= 1 V> V 

(D + 1)™ ^ ^ 

^ ' r— L q—r 

by equation (84). Second, one checks that P™ = for n < q 



E 

0<r<q 



D + l 

<7 



(15) 



(-1) 



q—r 



E (-l)^ r (^)^^ [xd x } n (l-x)"\ x=1 =0, for all n< g 

0<r<q 



(16) 



One thus has: 



(IrJ) = VnPS — = Vn 7 — -- V (-1 

V zw ^ D D + 1 ^ D + l)" ^ V 

n>0 V ; 0<r<D 



D-r 



\D-r 



D 



E (-D 

0<r<D 

(D+i) E (- 1 



E< 



D + l 



(17) 



D-r 



0<r<L> 



D 

r y (D + 1 - r) 2 



□ 

In particular for D = 2, 3, 4 one gets ( 1 7^ | ) = 9/2, 22/3, 125/12 respectively. 

Now let us make two remarks. Firstly, Lemma 2 declares that, on average, the depth of a vertex in a 
melonic D-ball is, up to a constant rescaling by Aa, just the tree depth in the associated (D + l)-ary tree. 
Secondly, from the iid tensor model, one has that the Feynman weight of the melonic two point graphs 
is equal to 1. Thus, the family of melonic D-balls corresponds to uniformly distributed trees. These two 
criteria set up an application of the non-trivial results established in [48] , which leads one to the following 
conclusion. 
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Theorem 1. Under the uniform distribution, the family of melonic D-balls converges in the Gromov- 
Hausdorff topology on compact metric spaces to the continuum random tree: 



mi 



A Ai /i5±D« 



(T2e,d 2e ). (18) 



D 



While we relinquish most details to Appendix C, it might be beneficial to explain at least the concepts 
involved in the statement of Theorem 1. 

One is familiar at this stage with the family of melonic D-balls. For the purposes of the theorem above, 
a melonic D-ball with n internal vertices is represented as a metric space: (m„, d mn / (Aa\/(-D + l)n/D)). 
Then, one looks at sequences of melonic -D-balls, with increasing number of vertices, such that the element 
of the sequence at any given n is chosen randomly with respect to the uniform distribution at that n. The 
theorem states that, in the Gromov-Hausdorff topology on these associated metric spaces, the elements of 
such a sequence converge to the metric space known as the continuum random tree [53]: (72e,^2e)- 

A melonic D— ball as a metric space: One should have a certain familiarity at this stage with the 
family of melonic D-balls. With this in mind, one denotes a random melonic D-ball with n (internal) 
vertices by M n . Then, one denotes this set of n vertices by m n . As always, these n vertices are in 
correspondence with the n (elementary) vertices of the associated rooted colored (D + l)-ary tree. 
Through this correspondence, there is a word associated to each element of m n . As a result, one can 
put a lexicographical order on the elements m n , that is, the vertices of m n are ordered according to 
how their associated words occur in the dictionary. With respect to this order, one denotes the rth 
vertex in m n by r for r £ {0, . . . , n — 1} (and its associated word by w(r)). Then d m „(fi, r 2) is the 
graph distance between r\ and r 2 in M n . Just to be clear, d„ ln (0,r) is the depth A(w(r)) defined 
earlier. 

The distance between any two vertices can be well estimated from the depth. Consider two vertices 
r\ and r 2 with words wuqu and wvqv. Thus, the two words have w in common, uq ^ vq are the 
first letters at which the two words differ and u and v denote the remaining letters in the words 
associated to r\ and r 2 . One denotes the vertices corresponding to w , wuq and wvq by r s , rf and 
rf, respectively. One has A(uqu) = d mn (rf, n) and K(vqV) — d mn (rf, r 2 ). 

The crucial point is that all the descendants of rf (including n) are connected to the rest of the 
melonic ball by a path going necessarily through one of the vertices in the distance array of rf 
(possibly r s itself). One calls these vertices r^' 1 , for i £ {0, . . . D}. The same holds for rf and r 2 , and 
one denotes by r^ 3 , j £ {0, . . . D} the vertices in the distance array of rf. The distance between rf 
and r*' 4 is at most 1. By the triangle inequality, in the triangle formed by r\, rf and rj' 1 , one has: 

d mn (ri ) rf)-l<d mn (n,r 1 ;i )<d mn (r 1) rj) + l Vi , (19) 

and similarly for the triangle formed by the three vertices r 2 , r| and r 2 J . The geodesic path from 
T\ to r 2 passes through some fixed r*'* and r s ^ 3 . Hence, for some fixed i and j, one has: 

d mn {ri,r s { % ) + d m „(r 2 ,r 2 u ) <d m „(ri,r 2 ) . (20) 

On the other hand, the path n — > rj 1 ' 4 —¥ rf — > r s — > r| —¥ r 2 J — > r 2 connects r\ and r 2 , hence 

d mn {ri,r 2 ) < d mn (r 1 ,r s 1 ' 1 ) + d mn (r 2 , r 2 a ) + 4 , (21) 

and one concludes that: 

dm n (ri,r{) + d m „(r 2 ,r|) - 2 < d m „(ri,r 2 ) < d r „„(ri,rf) + d m „(r 2 ,r|) + 6 

=> \dm„(ri,r 2 ) - A(uqu) ~ A(u i;)| < 6 . 



(22) 



To make M n a compact metric space, one needs a continuous metric. Thus, one must interpolate 
between the integer points on the integer grid (ri,r 2 ), for r 1; r 2 £ {0, ...,n — 1}. A piecewise 
linear interpolation on the triangles with integer co-ordinates (ri,r 2 ), (r± + l,r 2 ), (ri,r 2 + 1) and 
(n + 1, r 2 + 1), (n + 1, r 2 ), (ri, r 2 + 1) suffices. 

As n gets large, one would like to ensure convergence to some compact metric space (rather than 
just letting the structure get infinitely large). This requires a n-dependent rescaling of metric. This 
is the genesis of the factor Aa\/(D + l)n/D. Then, (m n , d mn /(A&yJ(D + \)n/D)) represents the 
melonic D-ball as a compact metric space. 
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A melonic D— ball as a random variable: It is worth noting that in the previous description, it was 
slipped in that M n denotes a random melonic D-ball with n vertices. The set up of the (tensor) 
model ensures that at every n, the set of melonic D-balls is endowed with a uniform distribution. 
Thus, one draws a random melonic Z?-ball from this set according to this distribution. Moreover, 
this entails that in this context, convergence means stochastic convergence, that is, convergence in 
distribution. 

Continuum Random Tree: A continuum random tree (CRT) (72e,d2e) is defined as a rooted real tree 
encoded by twice a normalized Brownian excursion e and endowed with a metric di e . 

One is probably more familiar with rooted discrete trees, of which the colored rooted (D + l)-ary 
trees are examples. Colored rooted (D + l)-ary trees (like all discrete trees) have an associated 
contour walk. Consider such a tree with n (elementary) vertices. (For simplicity, we shall consider 
its defoliated version, that is, all leaves removed). Starting from the base of the tree, one traverses 
the perimeter of the tree, passing from one vertex to the next in unit time-steps. One considers 
the following continuous function /(i), with /(0) = 0. As one travels, f(i) = d T (v) + 1, where 
v is the vertex one encounters at the ith time-step. (For the value at intermediate times, one 
linearly interpolates between the time-steps.) The procedure is illustrated in Figure 8. Given the 
construction, one has that the journey ends at time-step 2n, with f(2n) = and f(t) > for 
< t < 2n. One has thus associated to any tree some (fixed) walk /. For random trees with 2n 
vertices, the contour walk becomes a random walk with 2n steps. 





Figure 8: A defoliated (D + l)-ary tree and its associated contour walk. 

Any real continuous function f(t), such that /(0) = /(l) = and f(t) > for < t < 1, encodes a 
rooted real tree 7/. To get to the tree, one must set up the following equivalence. For all s,t € [0, 1], 

set m f (s,t) = inf max(Sit) < r < m i n(Sit) f(r). Then: 

s~t <i=> f( s )=f(t)=m f (s,t). (23) 
Then, the rooted real tree is the quotient: 7} = [0, 1]/ ~. The distance on the tree is given by: 

d f (s,t) = f(s) + f(t)-2m f (s,t). (24) 

One can pick out the branching vertices of the tree as those values in [0, 1] that are congruent to 
two or more other values. This real tree differs from a discrete tree in that one has precise distance 
information along the edges of the tree. 

The Wiener process is a stochastic process Wt (that is a random variable for every time t) such 
that Wo — 0, t — > W t is almost surely continuous, W t has independent increments and W t — W s is 
distributed on a normal distribution of mean and variance a 2 = t — s for s < t. The normalized 
Brownian excursion e t is a Wiener process conditioned to be positive for < t < 1 and be at at 
time 1. It is formally represented by a the path integral measure 



dfi e = ^ [dq(t)] 



e-^oBWl 2 *, (25) 

?(0)=?(1)=0 
q(t)>0 



with Z a normalization constant. 

The CRT (72 e , d^e) is the random tree associated to twice a normalized Brownian excursion 2e. 

Gromov-Hausdorff topology and convergence: Since one considers a sequence of random metric 
spaces, one should accurately define the space of metric spaces along with an appropriate topol- 
ogy. This is provided by the Gromov-Hausdorff topology on the space of isometry classes of compact 
metric spaces. 
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To begin, one considers a metric space (E,cIe). The Hausdorff distance between two compact sets, 
K\ and K 2 , in E is: 

dHauW^i, K 2 ) = M{r\Ki C K r 2 ,K 2 c #[} , (26) 
where K\ = {J xeK . BE{x,r) is the union of open balls of radius r centered on the points of Ki. 
Now, given two compact metric spaces (Ei,di), the Gromov-Hausdorff distance between them is: 

d GH (E 1 ,E 2 ) = M{d Haus{E) (MEi),ME 2 ))} , (27) 

where the infimum is taken on all metric spaces E and all isometric embeddings </>i and <p 2 from 
{E x ,di) and (E 2 ,d 2 ) into (E,d E ), 

It emerges that K, the set of all isometry classes of compact metric spaces, endowed with the Gromov- 
Hausdorff distance dan is a complete metric space in its own right. Therefore, one may study the 
convergence (in distribution) of K-valued random variables. 

Of course, the Gromov-Haudorff topology is not the exclusive topology for these metric spaces, but 
fortuitously, it is well-adapted to the study of quantities that are dependent on the size of the melonic 
D-balls, quantities such as the diameter, the depth, the distance between two random points and so 
forth. 

Superficially, the Gromov-Hausdorff topology appears to be quite an eyeful. However, convergence 
in the Gromov-Hausdorff topology is a consequence of any convergence of (some sequence of) E\ to 
E 2 embedded within a some common metric space E. In this case, the place of E\ is taken by the 
sequence of random metric spaces M n . Additionally, E 2 is the continuum random tree T 2e given in 
[53]. Then, the convergence claimed above rests on two points: i) Skorohod's representation theorem 
states that there exists a metric space fi, within which M n (for all n) and 72e can be embedded, and 
is such that the image of d mn approaches the image of d 2e almost surely as n — > 00; it) under the 
uniform distribution: 

( d mn ( Sl n,s 2 n) \ , « 

[ wra/sj (M)6M „ ™ <*•<-*»<«->*-.» ■ (28) 

This second point is the result proven in the appendix of [48] and to which we also devote Appendix 
C for explanation. It is clear that the involvement of Aa in the rescaling of the metric d mn is a 
highly subtle point. 

Importantly, the Gromov-Hausdorff distance between M n and T 2e is bounded from above by: 

sup { (d nn ( Sl n, s 2 n)/(A A ^(D + l)n/D) - d 2e (si, s 2 )) : (s u s 2 ) £ [0, 1] x2 } . (29) 
Thus, using (28), the stated convergence is ensured. 

With the previous explanation, it is now clear that the Hausdorff dimension of the family of melonic D- 
balls may be read off from (18) as the inverse of the exponent of n in the rescaling of the metric, that 
is: 

d H = 2 . (30) 



5 Spectral Dimension 

One uncovers the spectral dimension of a structure by analyzing an appropriate diffusion process on the 
structure in question. As a simple example, consider a diffusion process on a flat D-dimensional Euclidean 
geometry. One finds that the return probability attached to this process is: P{a) = a~ D / 2 7 where a is the 
diffusion time. One may extract the spectral dimension by taking the logarithmic derivative: 

alogcr 

Obviously, the spectral dimension coincides with the Hausdorff dimension in this elementary case, but this 
is not true in general. 

In the case of the branched polymers phase arising in Dynamical Triangulations (where one deals with 
uncolored binary trees). The appropriate diffusion process generates an average return probability, which 
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in turn gives rise to a spectral dimension ds = 4/3. This was shown by Jonsson and Whcater in [47], using 
an argument that we shall follow rather closely. 

It is worth mentioning, however, that for graphical structures, the spectral dimension depends rather 
strongly on the connectivity of the graph. This comes into play in the analysis below in that, a priori, 
one has a choice of graphical structure upon which one can place the diffusion process: the rooted melonic 
graphs, the melonic D-balls, the rooted colored (D + l)-ary trees or even the stack D-spheres. We 
shall choose the rooted melonic graphs over and above the others. Although the direct analogue of the 
binary trees are the (D + l)-ary trees, the rooted melonic graphs are the topological dual to the melonic 
D-balls and so more directly capture the connectivity of the manifold. Meanwhile, a diffusion process 
on the melonic D-balls themselves is difficult to analyze. For rooted melonic graphs, the appropriate 
diffusion process generates a return/transit probability, where transit refers to the process of traversing 
from one external vertex to the other. 2 The spectral dimension is, however, still extracted from the return 
probability. 

To proceed, consider a rooted melonic graph, with external color 0, contributing to the connected 2-point 
function. Such a graph is drawn in Figure 9. 

M 1 




M D 



Figure 9: A rooted melonic graph M. with sub-melons M. 1 . 

One denotes it by M.. Due to its iterative structure, such a melon is constructed from D + 1 rooted melonic 
graphs, each with a distinct external color, connected in the fashion illustrated in Figure 9. One denotes 
this property by: M = M 1 U M 2 U . . . M D U A4°, where M l labelled the rooted melonic graph with 
external edges of color i. 

Any connected 2-point graph, hence any rooted melonic graph, has two external vertices, one white and 
one black. Unless distinguished by a special name, the white external vertex is referred to as o, while the 
black one is referred to as •. For M. itself, they are called them respectively: o = I and • = O. 

First return/first transit probabilities: Consider a random walk on the melon M.. If the walker is 
at one of the external points, o = J or • = O, one can see from Figure 9 that it steps with probability one 
to its unique neighbor. Meanwhile, if the walker is at any of the internal (D + l)-valent vertices, it steps 
with probability j^ri to one of its D + 1 neighbors. 

One denotes by P)^(t) the 2x2 matrix encoding the first-return/first-transit probabilities, in time i, for 
the rooted melonic graph M.. In detail, its four elements are: 

- P]^° (t) = Pj^ 1 (t), the probability that the walkers starts from the external point o = I and returns 
for the first time to o = I at time t without touching the external point • = O in the intervening 
time. 

- P]^' '(t) = Pj^°(t), the probability that the walkers starts from the external point o = I and reaches 
for the first time • — O at time t, without touching / a second time. 

- Pm° {t) — P]^ 1 (t), the probability that the walkers starts from the external point • — O and reaches 
for the first time at o = I at time t without touching O a second time. 

2 To be more precise, a diffusion process is characterized by the diffusion equation along with appropriate boundary con- 
ditions. For the rooted melonic graphs, it is appropriate to choose cyclic boundary conditions. This allows the interpretation 
of the process as occurring on a closed manifold, that is, the related closed melonic graphs. 
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- Pjtf°(t) = P]^°(t), the probability that the walkers starts from the external point • = O and 
returns for the first time to • = O at time t without touching the external point o = 7 in the 
intervening time. 

Note that for any M, Pjj, M (0) = P^°*(0) = Pjfr°(0) = P^"(0) = 0. Moreover, the simplest melonic 
two point graph, denoted M(q), consists of exactly one line connecting the two external points. Its first- 
return/first-transit probability matrix is: 



Return/transit probability: Importantly the first-return/first-transit probability matrix of Ai de- 
termines the return/transit probability matrix of A4 in time t, denoted by Pj^(t). It is defined similarly 
to Pj^(t), except that the walker is allowed any trajectory between the end and final points. Indeed, a 
generic path from o = 7 to o = 7 can be decomposed as a word on I and O starting and ending with 7 
(for example, 77, 777, IOI and so on). Each pair of consecutive letters represent a walk from the first 
to the second letter that does not touch either external point in the intervening time. Similarly, a walk 
from o = 7 to • = O can be decomposed as a word on 7 and O starting with I and ending with O. Let 
us denote by w q , such words over 7 and O of length q, and denote w q (i) the ith letter of w q . Then, the 
return/transit probability in time t is: 

pxv(t) = 6 XY 6 t ,o + P)j i XY (t) 

(33) 

+ 2^2^ 2^ M {tol^M \PV---r M (t q -l)P M (t q ) . 

q=l W q toH \~tq=t 

At this stage, one introduces generating functions for both quantities, P]^ Y (y) — J2tV t PM <Y "(*) an d 
Pffi {y) = J2t y t P > M r (^)- The above relation then becomes: 

p^ Y (y) = s XY + P^iy) + EE PM Wq{1) (M)PM q(1)Wq(2) (y) ■ ■ ■ P^ q{q - 1)Wq{q \y)P^ q{q)Y iv) ■ (34) 

q=l w q 

Noting that w q (l) € {1,0}, w q (2) G {1,0} and so on, the above equation may be rewritten in matrix 
form (with obvious notation): 

p m (v) = i + Piiiy) + [Piiiy)} 2 + ■■■= 1 _p L(y) . ( 35 ) 

or in detail: 

Pm(v) P°M(y) 
P'Aiy) P'Aiy) 

l A - Pm'(v) Pm'Iv) 

[i - PTiv)} [i - Pjtrm - p 1 A'°(y)Plk°'(y) { pft°{ y ) 1 - p^°(y) 




Decomposition and iteration: The first-rcturn/first-transit probability matrix of Ai can be computed 
in terms of the first-return/first-transit probability matrices of its sub- melons M 1 , . . . , Ai D and Ai° . 

A generic walk contributing to (1,11) (that is, a walk starting at / and returning for the first time to I 
at time t, which does not touch O in the intervening time) can be further decomposed as a word over 
/, A, B, O (see Figure 9). Starting with I, the walker first steps to A. Once at A, the possibilities increase. 
The walker may travel out from A and return without hitting B, or may traverse to B. Moreover, if the 
walker is at B at some time t, she may travel out (also in the direction of O, although never hit it) and 
return, or may traverse back to A. These sub-paths may occur any number of times before eventually the 
walker returns to A a last time and steps back to /. Similarly, a walk contributing to (1,10) starts from 
/, passes through A, then A and B a number of times, before finally going to B a last time and jumping 
to O. Thus: 

(1,77) = 7 A...AB...BA. 
(1,10) = I A...AB...BA. 
(1,01) = O B ...BA...A I 
(1,00) = OO or O B...B A 



.AI 

AB...BO 



.AB...BO 



(37) 



1G 



While (1, OJ) is described similarly to (1,70), the are more walks than one might expect contributing to 
(1,00). Indeed, there arc walks, signified above by OO, that start from O, go into the melon M° and 
return to O without ever touching B. 

The first-return/first-transit probabilities between A and B may be written in terms of the first-return/ 
first-transit probabilities for the sub-melons: M 1 ,...,M D and M°. Indeed, for a random walk to go 
from say A to A it needs to chose one of the sub-melons M 1 , . . . M D and go from • to • in it. Thus: 

, (38) 
P MBA (t) = ^ 1 (P 1 M°'(t) + --- + P 1 M°'(t)) 

p hBB (t) = ^(pl 4 r(t) + ---+p 1 M o °(t)+p l M o °(t)) , 

where again one notes that the first return probability for B to B is special, as the walk may go also into 
the melon M°. (Note that on the right hand side of equation (38), • = A and o = B.) Denoting by P 1 (t) 
the first-return/first-transit probability matrix between A and B, one has: 

PH*) = nTlPUt) + -.. + j+r^t) + ^ J) . (39) 

Furthermore, the first-return/first-transit probabilities between B and O are: 

P hOB (t) = Pfc°(t) , P^°(t) = ^P^o'it) , P^ 00 (t) = P^ '(t) . (40) 

On this occasion, let us denote by w q the words of q letters over A and B. The walks of first-return 
starting and ending with I decompose as: 

+ - j"EE E ^P 1 ^^( 1 )(i )P 1 ^( 1 )^( 2 )(t 1 )...P 1 ^ 9 ('?-l)^(9)(i 9 _ 1 )pl^(9)^(t (J )^ . 

q=l Wq t H \-tq=t-2 

The first step of the walk is from I to A. At the second step, the walker either returns with probability 
(D + 1) _1 to I (the first term) or proceeds down one of the sub-melons M 1 , . . . , M D . The walker can 
then either return to A without touching B in t — 2 steps and go to I with probability (D + l) -1 at the 
last step (the second term), or go from A to B a number of times, end in A and jump with probability 
(D + 1) _1 at the last step in I (all the other terms). Similar considerations lead to the following equations: 

= E p v ' AB (to)p v ' BO (h) 

t +ti=t-l 



+EE E (p 1 ^^^)..^ 1 ^^ 5 ^)^ 1 ^ ^!)) , 

q=l w q t H K„+l=t-l 

P^if) = E ^ 1;OS (io)P 1;SA (ii) 

t +ti=«-l 

^ OO 

+ittiEE E (p liOB (<*.)P lifl -W(ti)...p li "«w A (^i)) , 



(41) 



q=l w q t H ht 9 +i=t— 1 



P#'(t) = P^°°(t)+ ^ P 1 ^^)^ 50 ^) + E ^ 1;OB (to)P 1;BiJ (il)P 1;BO (i2) 

t +ti=« i +ti+t2=t 

oo 

+EE E (p 1;OB (to)p 1;B ^ (1) (*i)---p 1; ^ (9)B (i 9+ i)p 1;SO (^ + 2)) • 

<j=l w q t H h*<,+2=* 
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As one might expect, P^ (t) requires special consideration. The first term, P 1; °°(i), represents the walks 
which start from O and end in O without touching B. As before, the equations simplify for generating 
functions: 



PmUA DTi y2 + DTi y2pl '' AA{y) 

1 CO 

y 2 ^ p^Aw q {\) (yjpl;w q (l)w q (2) ^ _ _ _ pl:w q (q)A^ 



+ 

~D 



9=1 



Ptf'iv) = yP 1 ' AB (y)P 1 ' BO 



+y^2 P 1 '^"^ {y)P 1 ' w ii 1 ) w ii 2 ) (y) . P 1 '^"^ 8 (y)P 1 ' BO (y) 
9=1 

' yP^ OB {y)P^ BA {y) 



D + l 



+ _i P VOB (y)P v < Bw «W (y)P 1 ' w "^ w " ( - 2 \y) . . . P 1 ' w "^ A (y) 



9=1 

Ptf'iv) = P 1;00 (!/)+P 1;0S (!,)P 1;fl0 W+P 1;0B WP liBS (# 1;B0 (!/) 

oo 

+ P 1;OB (y)P 1;B ™*( 1 >(y)P 1; ™<>W , ''<>( 2 )(2/) . P 1;M, <^ B (y)P 1;BO (y) 
9=1 

Furthermore, in matrix form, they look yet simpler: 







P^d/^lO P 1;OB (y) 



y o 



£>+l 





V.BOl 



pWV{y)) ' 



where: 



a = 



1 N 

1 0/ • 



(42) 



(43) 



(44) 



Substituting the various probabilities as a function of the sub-melons one gets a recursive equation along 
with an initial condition: 



Pit(y) 



o o 



1 (y 



P£o'{y))^D + l V0 P^ °(y) 



x a 



1 "d + iK i(!,) + -^ d(! ' ) + 



<o°W o 





"(o P^y)) ' (45) 



y 
V 



In principle, this solves the problem of determining the first-return/first-transit probabilities for an arbi- 
trary melon. For example, the first non-trivial melon is the elementary melon of color 0, denoted by A4(i), 
and one gets: 



Pm (1) (v) 



1 



y 
D + 1 V° V 



1 - 



D 



y 



y o 

y 



-=— v 2 



( 1 



1 - 



( D +l)2 2/ 2 Vd + i^ 



p + i 

Defining some auxiliary 2x2 matrices cleans up the formulae quite significantly: 

E% = 5 a a 5l where a, b, a, e {1, 2} . 



Thus: 



P 1 



= P 22 P_t, P 22 + (p 12 y + P 22 P^ P n ) 



P 1 



y 

y o 



ya- 



rn 



(47) 



(yP^+P^oP 22 ) (48 ) 
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Note that by induction one trivially obtains that the first-return/first-transit matrix is symmetric 



p 1 M'°(y) = P 1 M'(y)- (49) 

Ultimately, one must solve this system in order to obtain a closed form for the return/transit probability 
and thereafter extract the spectral dimension. Suggestively, the matrix form of the system in (46) has 
superficial similarities to the system given in [47]. However, in comparison to [47], there are two added 
complications: i) it is a matrix rather than a scalar equation; ii) the final term in the denominator 
obscures a direct recursive solution. To circumvent these problems, one considers first so-called simple 
melons. Later, we shall argue that the general case is just a small correction that does not affect the 
spectral dimension. 



Simple Melons: Simple melons are those rooted melonic graphs M. conditioned by: i) Ai° = and ii) 
for M any sub-melon of color i in A4, then M l = 0. In terms of the associated (D + l)-ary tree, a simple 
melon corresponds to a tree such that no branch possesses two consecutive lines of the same color. Then, 
the final term in the denominator vanishes and the formula reduces to: 



i 



(50) 



M ' 



Lemma 3. Pj^ = a + bo for all simple melons A4, where a, b € R and a implicitly multiplies the 2x2 
identity matrix. 

Proof: First of all, PL = ya, so it is of the claimed form. Then, one uses an inductive argument. One 



(0) 

notices the following general matrix relationship: 

a 



a(a + f3o) 1 o = er(- 



P 



.. P 2 a 2 -/3 2 ' a 2 -P 2 a 2 ~p 2 ' 

Thus, should each PL ; be of the claimed form, one can use (51) to rewrite P)^ in the same form. 
Noting moreover that o~(a + ha)a = a + ba, it follows that the recursion may be rewritten as: 

1 



(51) 



□ 



and that all the PL; 



bid may be diagonalized simultaneously in the basis: 



(52) 



with eigenvalues X\ 



M*;(l,2) 



= a, ± bi, that is: 



V2V 



1 

71 



p 1 



1 1 1 



2 VI -1 



\1 

A M;1 




\1 

A M:2 



1 1 

1 -1 



(53) 



In its diagonalized form, the recursion is: 



\i 

A M;1 



(. 



y 



A M:2 



D + l-Ei Al Mi;1 




V 



D 



A M (0 y,l = y> A M {0 y,2 = ■ (54) 



Thus, one may restate the scalar recursion for each eigenvalue in terms of a single \m{u)'- 



A M (y) 



D + i-Ei A M<y) 



and X M(0) (y)=y 



(55) 



with Aj^.j = \m(u) and = a m{~ y)- At this stage, the system for Xm(v) has almost coincided with 
that of [47]. However, one must still tread carefully: first the powers of y differ with respect to [47] and 
second, one must deal with the extension to D + 1 colours. The argument is checked below. First, one 
introduces the function: 



hM(y) 



-(l-A M (y)) , 



(56) 
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so that the recursion may be rewritten once again as: 

i + y + EihM*(y) 



h M(y) 



and h 



M (n) (y) = j— (l-y) 



i + (i-y)Ei^*(w) 

One defines the generating function: 



(57) 



(58) 



One can show that Q(z, y) has a simple pole at y = 1 by showing that the sum J2m zP /^m(^) converges. 
At y — 1, the recursion and initial condition are given by: 



h M (X) = 2 + ^2h M i(l) and h M(0) (l) = l 



(59) 



If the melon A4 has 2p internal vertices, the contribution of the pole to the sum is: 

h M (l) = 2 + ^[(Z> + l) Pi + l] = 2 + D + (D + l)J2Pi = (D + l)p + l, 



(60) 



where one takes into account that the melon M. with 2p vertices satisfies p = J^iPi + 1, if the sub-melons 
M. 1 have 2pi vertices. Thus, 



^-y)Q{z,y)\ v=1 =^C p \ where C p = ^ (° P '+ 1 

p— 



(61) 



which converges for z < 



(D-l 



,D-1 



Z . 



The information about the spectral dimension is contained the remaining non-pole part, so one defines: 

d 



Q(z,y) = -~(l - y)Q(z,y) - E^/^^S ~h M (y) 

Rather than tackle Q(z,y) all in one go, one concentrates for a moment on Q n {z), where: 

L — ^ TJ I 



(62) 



(63) 



so that: 

Defining further: 
one can can express: 

n+l 

QnW = E(- 1 ) , ' +1 
r=l 



dy" 



~ d n+ 1 



M 



h M = d^i h M{y) 



E 



oi,...,a„+i 



...[(n + l)!] a "+iai!...a„+i! 



E 

M 



■AO) 

1 M 



A3) 
r+1 



(64) 



(65) 



(66) 



where more details are given in equations (47) to (52) of [47]. From this point on, we change notation 
slightly. We shall use Nm to denote half the number of vertices in A4 , while p will represent the number of 
indices in subsequent formulae (we do this to keep notational similarity to the analogous argument given 
in Section 5 of [47]). 

We denote: 



(67) 



M 
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In order to derive the asymptotic behavior of Q n (z) we first note that 



M \ ll M) / M ' \ a M) 

= Y, zNm Y[( h Mr = i3-( ai ® 1 -- a "+ 1 ® n+1 > , CH®i = i,...i (68) 

where we used (60) and (65). Thus the asymptotic behavior of the term of order r in (66) can be obtained 
by integrating r + 1 times the asymptotic behavior of 7j( a i® 1 > - >°n+i® rl + 1 ). We show in Appendix D that: 

Lemma 4. For all ni, . . . n p i/ie following asymptotic behavior holds 

H {n u ...n P )^ ~^Z Nm h$ ) h ( tf ) ...hfc ) ~ U 3-P-|(»i+n2+-+»*) ; fls 2 j Zq _ ( 69 ) 
M 

with u = 1 — z/zq . In particular 

jj(ai®l,...,a n+1 ®n+l) ^ u |-r-f(n+l) ^jq-j 

Substituting the leading order behavior (70) and integrating r + 1 times, one obtains the leading order 
behavior: 

Q {z) ~ log(l - z/zq) , Q B («)~(l-z/zb) _|n , Vn>0. (71) 

Consequently, the most singular part of Q(z, y) as z f z is the sum of logarithmic piece and a function of 
(1 — y)(l — z/zq)~z. Thus, one arrives at an expression which is identical to equation (57) in [47]: 

oz l~z/z \(l-z/zo) 5 / 

for some function $. The rest of the analysis coincides exactly with the one given in [47] section 4.2, up 
to the two comments made below, so that the result for the spectral dimension is: 

ds = \ ■ (73) 



There are two comments to make here. The first point to note is that the dQ/dz obtained in (72) pertains 
to the system for Ax (y) with initial condition Xm (0) (v) — V-> that is, to A^ ;1 (y) . One should recall that this 
scalar problem arose via the diagonalization of a matrix problem. So, a priori, it does not simply provide 
the behavior of the return probability, but rather the sum of this and the return/transit probability. To 
remedy this, one should repeat the process for the system with initial condition Awi (0) (y) = —y. The 
resulting Q(z,y) has a pole at y = — 1 and after removing this, one finds that the same procedure leads 
to (72), except with y — > —y on the right hand side. This provides the behavior of the difference of the 
return probability and the return/transit probability. By following an argument analogous to that given 
in Section 4.2 of [47], for the quantity: 

oz dz 1-z/zo \(l-z/z a )2/ 1-z/zq \(1- z/zq)* J 

in the regime \y\ < 1, one extracts the behavior of the return probability, and the spectral dimension stated 
in (73). 

The second comment refers to the non-simple melons, which have been neglected up to this point. We 
argue rather indirectly that their inclusion has a negligible effect on the spectral dimension. Durhuus, 
Jonsson and Wheater have shown in [55] that the spectral dimension of the infinite Galton-Watson tree 
is 4/3. We explain in Appendix C that rooted melonic graphs correspond to such objects in the infinite 
limit. To establish ds = 4/3 using this line of reasoning, one would have to have to repeat that analysis 
again with color information and so forth. Although we do not do this in detail, the results we provide 
here overwhelmingly back up this claim. 
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A The road to melonic graphs 

Let us construct the class of independent identically distributed (iid) models. We shall attempt to be 
precise without being especially detailed. We refer the reader to [56] for a more thorough explanation. The 
fundamental variable is a complex rank-D tensor, which may be viewed as a map T : H\ X • • • x Hp — > C, 
where the Hi are complex vector spaces of dimension 2Vj. It is a tensor, so it transforms covariantly 
under a change of basis of each vector space independently. Its complex-conjugate T is its contravariant 
counterpart. 

One refers to their components in a given basis by T n and T n , where n — {n 1; . . . , njj}, n — {n\, ...,«.£>} 
and the bar (— ) distinguishes contravariant from covariant indices. 

As one might imagine, with these two ingredients, one can build objects that are invariant under changes 
of bases. These so-called trace invariants are a subset of (T, T)-dependent monomials that are built by 
pairwise contracting covariant and contravariant indices until all indices are saturated. It emerges readily 
that the pattern of contractions for a given trace invariant is associated to a unique closed D-colored 
graph, in the sense that given such a graph, one can reconstruct the corresponding trace invariant and 
vice versa. 3 

In Figure 10, we illustrate the graph B\ : the unique closed D-colored graph with two vertices, which 
represents the unique quadratic trace invariant trgj (T, T) = T n S nfi T Tl . 



tl£>i (^3 ^0 T ni j l27l3 S nirll S n2rl2 fin-jris T niri2ri3 

Figure 10: A closed Z?-colorcd graph and its associated trace invariant (for D = 3). 

More generally, we denote the trace invariant corresponding to the graph B by trg (T, T) . 

From now on, we shall make two restrictions: i) all the vector spaces have the same dimension TV and ii) 
we consider only connected trace invariants, that is, trace invariants corresponding to graphs with just a 
single connected component. 

Given these provisos, the most general invariant action for such tensors is: 

oo 

where is the set of connected closed D-colored graphs with 2k vertices, {t^} is the set of coupling 
constants and oj{B) > is the degree of B (see [1] for its definition and properties). This defines the IID 
class of models. 

3 While we refer the reader to [1] for various definitions, it is perhaps not unwise to match up right here the defining 
properties of a closed D-colored graph B with those of a trace invariant: 

- B has two types of vertex, labelled black and white, that represent the two types of tensor, T and T, respectively. 

- Both types of vertex are D-valent, with matches the property that both types of tensor have D indices. 

- B is bipartite, meaning that black vertices are directly joined only to white vertices and vice versa. This is in 
correspondence with the fact that indices are contracted in covariant-contravariant pairs. 

- Every edge of B is colored by a single element of {1, . . . , D}, such that the D edges emanating from any given vertex 
possess distinct colors. This represents the fact that the indices index distinct vector spaces so that a covariant index 
in the ith position must be contracted with some contravariant index in the ith position. 

- B is closed, representing that every index is contracted. 
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A.l Free energy and melonic graphs 



The central objects for further investigation emerge from the free energy (per degree of freedom) asso- 
ciated to these models: 

£({M) = -^log(/ dTdf e"^- 1 ^)) . (76) 

When facing such a quantity, the standard procedure is to expand it in a Taylor series with respect to the 
coupling constants and to evaluate the resulting Gaussian integrals in terms of Wick contractions. It 
transpires that the Feynman graphs Q contributing to E({t&}) are none other than connected closed 
(D + l)-colored graphs with weight: 

where sym((?) is a symmetry factor, Bi p \ runs over the subgraphs with colors {1, . . . D} of Q and w(G) is 
the degree of Q. A few more words of explanation are most definitely in order here. The graphs Q € r( D+1 ) 
arise in the following manner. One knows that a given term in the Taylor expansion is a product of trace 
invariants upon which one performs Wick contractions. For such a term, one indexes these trace invariants 
by p € {1,... , p max }, that is, we index their associated D-colored graphs B< p )- A single Wick contraction 
pairs a tensor T, lying somewhere in the product, with a tensor T lying somewhere else. One represents 
such a contraction by joining the black vertex representing T to the white vertex representing T with a 
line of color 0. Thus, a complete set of Wick contractions results in a closed connected (as one is dealing 
with the free energy) (D + l)-colored graph. A particular Wick contraction is drawn in Figure 11. 




Figure 11: A Wick contraction of two trace invariants (for D = 3). The contraction of each (T,T) pair is 
represented by a dashed line of color 0. 

It requires a bit more work to reconstruct the amplitude explicitly, see [56]. Importantly, oj(Q) is a non- 
negative integer and so one can order the terms in the Taylor expansion of (76) according to their power 
of 1/N. Quite evidently, therefore, one has a 1/N— expansion. 

In the large-iV limit only one subclass of graphs survives, containing those graphs with lu(Q) = 0. In [27], 
it was shown that the only graphs of degree zero are the melonic graphs. Moreover, at leading order in 
1/N the graphs contributing to the two point function are the rooted melonic graphs. 

B Proof of Lemma 1 

Lemma. One has: 
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Proof: Note that: 



m +••• n q — n 



q = 



E 



q / «. ni + ...n r -7i 



rii!...n ! ^-^ \r J ' m!...n r ! * — ' \r 

jii *...n«>0 y r— 1 v 7 n 1 ,...n T .>l ' r— 1 



(78) 



Let us define L as a lower triangular q x q matrix with non-zero entries: L rs = rj for q > r > s > 1. 
Then, the S 1 ™ given in equation (7) provide the solution to the system of equations: 



r n = L rsS" , for all r such that 1 < r < q 

S 

We shall show that the lower triangular q x q matrix P, with entries: 



Pre = (-IT 



, for q > r > s > 1 



is the inverse of L. The key is to consider the following two series expansions: 



s=0 



(1 + x) 



1+ 



Taking their Cauchy product, one finds: 



(i + *r ^ = (i + *) r (H^+t = E * r ' 



E 

t<s<n 
t—s<r 



(79) 



(80) 



(81) 



,(-1) v . 
n — s I \t 



(82) 



However, from a direct expansion of the left hand side of (82), one can see that for t < r, the coefficient 
of x r in (1 + x) r ~ t ~ 1 x t is exactly 0. Hence: 



yt < r, 



E 

t<s<r 



1 V-iy-^l 







(83) 



As a result, one has the following two relationships: 



E L r> p * = E Q(- i r t (3 = | 



r < t 

1 for r = t 
r > t 



(84) 



Thus, one finds that as q x q matrices, L and P are inverse and the claim follows. 



(85) 



□ 



C Explanation of Theorem 1 in Section 3 

In Section 3, we provided the statement: 

Theorem. Under the uniform distribution, the family of melonic D -balls converges in the Gromov- 
Hausdorjf topology on compact metric spaces to the continuum random tree: 



A A ^/(D + l)n/D / 



"> (T 2 e,d 2 e) 



(86) 



As stated earlier, an exhaustive proof along the lines of [48] would inevitably be lengthy. Here, we shall 
give a very brief sketch. From the discussion at the end of Section 4, the proof of this result boiled down 
to showing that: 



d m „(sin, s 2 n) \ 



> {d 2e ( Sl ,s 2 )) {suS2)e[01]> 
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In fact, the proof of this point is not so direct. Rather, it rests on the following seminal result of Aldous 
[53]: 

Under the uniform distribution, the family of trees associated to a critical Galton-Watson process with 
variance a converges in the Gromov-Hausdorff topology on compact metric spaces to the continuum random 
tree: 

(t„,4=) >(T2e,d 2e ). (87) 



A Galton Watson process is a stochastic process {Zk}, evolving according to Zk+i = Y^f=i fj^ where 

(k) 

the Q are random variables, taking values in Z + , independently and identically distributed according 
some distribution [i. Colloquially, k is the generation number, Zk is the number of offspring in the fcth 
generation, while /i is the offspring distribution, so that p,j is the probability that any member has j 
children. A process is critical if its mean is one: ^Zj—aJ^j — 1- The variance of the distribution is 
a = Sj^o j(j ~ 1) Mj- Obviously, if Zq = 1, the process has one initiator and any instance of the process 
can be mapped to a tree. 

A Galton— Watson tree, T, is a tree generated by this process. For the theorem above, one is interested 
in Galton-Watson trees conditioned to have a total of n progeny, that is, trees T n with n vertices. Alterna- 
tively, critical Galton-Watson trees can be obtained as simply generated trees. The offspring distribution 
induces a distribution for T n on the set of all (plane) trees with n vertices, denoted T n . Specifically we 
associate to every tree a weight XV^{T n ) :— YivtT f-K? where v are the vertices of T n and k v the number 
of offspring of v. One then picks at random plane trees with probabilities proportional to this weight. 

P(T n ) = ^ * ,„ M T n ) . (88) 

The strategy is to show that the family of mclonic D-balls (seen as random metric spaces under the 
uniform distribution) correspond to some family of Galton-Watson trees (seen as random metric spaces 
under the P(T n ) distribution). 

The main flow of the argument has three parts. 

Part 1. From melonic D— balls to Galton-Watson trees: Fortunately, one knows already that mel- 
onic D-balls are in correspondence with colored rooted (D + l)-ary trees. The aim is to show that 
these are generated by a critical Galton-Watson process. The appropriate offspring distribution is: 
fiQ = D/(D + 1), (iD+i = l)j with the rest zero. Such a tree with a total of (D + \)n + 1 

vertices (n internal and Dn + 1 boundary) is denoted by Tm+iw+i- 

Part 2. Distributions: The weight IT M on such trees satisfies n M (T (D+1)n+ i) = D Dn+1 / (D + l) i - D+1 > l+1 . 
Thus, it is constant across T(D+i) n +\, the set of colored rooted (D + l)-ary trees with (D + l)n + I 
vertices, and corresponds to the uniform distribution: 

^<W + .) = ^, w„e re c;<-^ t7JT ^ TT (< D + ;>" + 1 ). m 

Cn D+lS> counts the number of elements in the set Tm+t)n+i- 

Part 3. Metrics: The problem becomes yet more nuanced when one moves to the metric spaces asso- 
ciated to the melonic D-balls. The metric space (m„, d mn /A^i l /(D + l)n/D) takes into account 
the n internal vertices of M n only. These are in correspondence with the internal vertices of some 
element of 7(D+i)n+i 5 sav ^Vu+iln+i- One denotes the defoliated tree corresponding to Tm + x\ n+ x by 
T n . One can cut through a lot of red tape by noticing that these defoliated trees are also generated 
by their own critical Galton-Watson process with binomial offspring distribution on the first D + 2 
weights: fij = ( D + l )D D+1 -i /(D + 1) £,+1 , for < j < D + 1 and the rest of the weights zero. This 
is critical with variance a = D/(D + 1). Thus the appropriate result to quote is that: 

Under the uniform distribution, the family of trees associated to a critical Galton-Watson process with 
variance a converges in the Gromov-Hausdorff topology on compact metric spaces to the continuum 
random tree: 

7f=k ,„ ) ► (T 2e ,d 2e ) . (90) 
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where dr n is the tree distance in the defoliated (D + l)-ary tree T n . The vertices of T n have a 
lexicographical order r € {0, . . . ,n— 1} generated from their associated words. Thus, for two vertices 
T\ and r 2 (in correspondence with two vertices ri, r 2 in the melonic -D-ball), rfx ii (ri,r 2 ) is the tree 
distance between r± and r 2 in the tree T n . 

One extends dx n to a continuous metric by interpolating between the integer points in the same 
fashion as for d mn . As a result of this (and so-called tightness of this family of rescaled tree metrics) , 
one has: 

/ dT n (sin,s 2 n) \ 



(si,s 2 )e[o,i] x2 



> (d2e(s 1 ,S 2 )) {si S2)e[01]x2 . (91) 



\y/(D + l)n/D 

Therefore, the final point is to show the following convergence of rescaled metrics: 

dm n ( L s i"J , L s 2«-J ) d Tn ( \_sin\ , ls 2 n\ ) 



A A ^/(D + l)n/D y/(D + l)n/D 



. (92) 



for all (si, s 2 ) e [0, l] x2 . (p) indicates that the results holds with probability close to 1. The route to 
this result is far from short. The (very) rough gist of the argument goes as follows. Consider some 
fixed pair (si, s 2 ) with s± < s 2 . Then, for each tree T n in the sequence, one can decompose the paths 
from the root vertex to u;(|_sinj) and w([s 2 nj) as: 

W([sin\) = W tru nk 'o 'branch and w(L«2«J) = Wt rU nk ^0 ^branch , (93) 

where U?trunk is the word corresponding to the common part of their ancestry, Iq, r$ are their respective 
first letters after they diverge and /branch, '"branch are their respective remaining letters. 4 

One now considers the following sequences of random variables: 

X n = ( Wtrunk, 'branch, ^branch, 'o, r o), and X n — (w trunk, 'branch, ^branch, 'o, r o) (94) 

where all the components are independent random variables in their own right and the tilde indicates 
that the variables are drawn conditionally on some values for the word-lengths: |u>trunk|i |'branch| 
and ^branch | at each n. The first three in each set are drawn from Wd+i, while the final pair are 
drawn from Id+i — {( a , b) : a, b € {0, ...,D} and a < b}. Furthermore, consider the following 
sequences of functions of these variables: 



<7n(A n ) — 1 / (■ r^_Li) n (I'branch | , |^branch| , A(^branch), A(rbranch)) 

~ - oir ~ (95) 

9n(.Xn) — \l TjJ+lJn (['branch | , |^branch| , A(/branch), A(rbranch)) ' 

For a moment, assume that one can show that as n — > oo, g n {X n ) — > {a Sl ,a S21 K^a Sl , K^a S2 ) for 
some a Sl and a S2 . Recalling eq.(22) and taking into account that |A(Zo'branch) — A(?branch)| < 1, the 
following inequalities hold: 

Mr„(|_Sl«J, L s 2"-J) - | 'branch | ~ Kranch|| < 2 

, (96) 

|rfm„(L s l n J, L s 2™|)| - A(/ b ranch) ~ A(r bran ch)| < 8 . 

(p) 



These imply that ^ (rj^ > i)„ "|AA^T„ (L^i^J , L S 2^J) — d mn ( \s\n\ , [s2«J)| — - — > and so the stated 
convergence (92) holds. 

There is, however, one final catch. One first shows only g n (X n ) — > (a Sl , a S2 , Aao Si , K.&.a S2 ) for 
some a si and a S2 ; rather than g n (X n ). By [54] (Lemma 16 therein), one implies the other if the 
distributions governing X n and X n converge as n — > oo. This element of the argument is proved in 
detail in [48] (Lemmas 35 and 36 therein). 



4 The notation is meant to be somewhat suggestive. If s\ < S2, the lexicographical order on the vertices ensures that 
u>([sinj) occurs to the left of w([s2n\) in the plane drawing of the tree. Hence, /branch and branch denote left branch and 
right branch, respectively. 
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D Proof of lemma 4 



In order to prove that: 

H (nu-n P )^ z) := J2z Nm h% l) h < tf ) ...h ( K ) ~ u I-p-|(»i+"»+-+"p) , as ztz , (97) 



M 



with u = 1 — z/z , one uses an inductive argument. 

To proceed, one requires a set of initial cases that are explicitly shown to satisfy the claim. Here, these 
turn out to be given by the set of #(°>— ,o). m other words, for any p, where rii = for all 1 < i < p. It 
emerges that: 



M 



= j2i(D + i)N+ir 



1 / DN + 1 



N 



DN + 1 V N 



E 

N 



N 



(98) 



since for a generic melon, M., with 2N vertices: h^} = (D + 1)N + 1. One has furthermore that: 



ho = Yz Nm =Y — 1 — ( 

^ ^ DN + 1 V 



X iV 

[JfO]* ~ constant + u 1 ^ 



ZW + 1 



z w ~ constant + h 1 / 2 



(99) 



#0 = 1 + *(#<))*> (l - £te[JT()p-i) 



[#()]' 



,.1/2 



To help with the inductive process, one expands ftW as a function of its sub-melons M 1 . The label r 
refers to the number of derivatives that act on the denominator, which yields: 



h 



M 



i + (i-y)Ei^«(i/) 



9"" 



l + y + J2h Mi (y) 



n n 

= 2*„, +<w+E^+ 2 E E rf « HyE^ 1 

i s=l ai...a„ j = i j 
T, a j= s ; T.j a i= n 

n—1 n—1 

+E E 4-1,.,.,) II [j'E* i; 

s=l ai...a„_i j = i j 

+EE^ r) E E *<r,.,. l'lyE'*..' 



r=l 



= 1 ai...a r 



where: 



d, 



r (-l) s+1 s! 



(100) 
(101) 



K rJ (l!) fl i ... [r\] a -a 1 l...a r \ 
and equation (66) has been used to obtain the final form. 

While the stage has been set to invoke the inductive hypothesis, it is beneficial to first tackle a simple 
example. Actually, this transpires that it is remarkable indicative of the general argument. The case in 
question occurs at p — 1, with m = 1. Evaluating the recursive equation (100) in this case, one finds: 

-1 2 



C-i + E^ + ^e 



5>S 



(102) 
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which leads to the following equation for ijW (recall that N = ^ Ni + 1): 

hV>(i-zD[hQ\ d - x ) 

= HO + 2zDH(°> [H^] - 1 + zD(D - l)H^H^ [H^] D - 2 + zDH^ [ff 0]D-i ( (103) 
with the use of the helpful remark that: 

Y, zNMh M' = zH^iHOf- 1 . (104) 

M 

The most singular term is the one involving H^ 0,0 \ thus one finds that the behavior of is: 

~ u- 2 , (105) 

which is in agreement with the claim. More importantly, one should note that the relation (103) allows 
one to express in terms of Hs with lower valued indices, although the number of indices may increase. 
Obviously, any inductive argument requires that one put an order on the set of .H"( ni >---> n 3>) . The above 
example should make one aware that this ordering is subtle. 

(As a brief aside, one notes that equation (103) is slightly different to the analogous equation given in [47] 
for the case of branched polymers. Obviously, appearances of D are to be expected. However, one may be 
slightly puzzled as to the presence of factors of z. There is a simple reason for this occurrence. In the case 
of binary branched polymers (trees) , the terms are labelled according to their total number of vertices R. 
This satisfies R = R\ + R2, where the Ri are the two sub-branches. In the case of rooted melonic graphs, 
the formula for N is: N = 1 + J2i N, as stated above equation (60). This extra 1 on the right hand side 
leads to the presence of the z— factors.) 

Returning to the main theme, the aim now is to put an appropriate order relation on the possible lists of 
indices (non- negative integers) S of H . Let us start by ordering the elements of S in decreasing order 
from left to right, that is S(i) > S(i + 1) if S(i) denotes the i'th element of S. Furthermore, one denotes 
by \S\ the number of entries in the list S. One orders the lists in lexicographical order, that is S a > Sb if 
one of the three following statements holds: 

. 5„(1) > S b (l). 

• S a (i) = S b (i) for 1 < i < i and S a (i ) > S b (i ). 

• S a (i) = S b (i) for 1 < i < \S b \ and |5 a | > |5 b |. 



Note that > is a total order relation: for any S a 7^ S b , one either has S a > S b or S b > S a . One states 
S a D S b if all the elements in the list S b together with their multiplicities are also elements in the list S a . 
If S a D S b and |5 a | > \S b \ then S a > S b . 

One uses this order relation in the following way. Consider an abridged form of equation (100): 

f& = Y^h^ t +T{n-l), (106) 

i 

where T(n — 1) contains the remaining terms, which by design possess derivatives of order at most n — 1. 
Let: S = {n, ...,71,711,... ,n ps }, with n > n± > ni ■ ■ ■ > n ps . One may write any H as: 



M 



His) = Y,z^[h<$ri[h 

M j 



Mi j i 

The generic form of a term on the right hand side of (107) evaluates to: 

Y[HM (108) 

i 

where Si denotes the set of indices pertaining to the sub- melon M 1 in the product. There are a number 
of possibilities: 
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- There is a term for each fixed i corresponding to choosing h^} in all the terms. It has Si = S and 

Sj = for all j 7^ i. These terms are brought over to the left hand side of equation (107). 

- The rest of the terms may be divided in several classes: 

• Terms with no factor T. They subdivide as: 

- no single set Si possesses all q indices of value n. Hence Si < S for all i. 

— a single set Si possesses all q indices of value n. But then Si must have less than ps elements 
drawn from {m, . . . , n ps }. One has Sj < S for j ^ i and S D Si, \S\ > \Si\ hence Si < S. 

• There is at least one factor of T(n — 1) and any number of factors of T(rij — 1) for the various 
j. In this case, any given Si contains fewer than q indices of value n, hence Si < S for all i. 

• There is at least one factor of some T(rij — 1) but no factors of T(n — 1). In this case, a given 
Si may contain the index n up to q times, at most ps — 1 indices {n il7 . . .n ik } drawn from 
{ni, . . . , n ps } with k < ps — 1 and the rest drawn from integers strictly smaller that 

max n , (109) 

ne{ni,...,n PS }\{jij 1 ,•■•«»,.} 

hence Si < S. 

One now proceeds to a more refined analysis of the leading divergent behavior of the various terms involved 
in eq. (107). One first associates a number, called the naive index do, to every hj$: 

d Q {h^):=-l-\n, (110) 

The definition readily extends to products, do^Ilj ^jw^) = (^m^ ■ ^ is convenient to regroup 

the terms in eq. (100) as 

i i r=l ^ ' 



where the terms in T' have either i) naive index at least — 1 — |n+|orn) naive index — 1 — |n + | but 

^ and f^j 



possesses at least two entries /iW and /iW corresponding to two distinct melons. It is time to proceed 



to the inductive argument. 

Inductive hypothesis: For all R — {n, . . . r t } < S the actual degree of divergence d a of (that is, 

HW ~ u d °) is: 

da(H^) = 1 -+d Q (f[h^). (112) 

Inductive step: Assume that the claim holds for all sets R < S = {ni, . . . , n p }. First, one expands H^: 

H (s ) = j2 zNm n(E* , +EE^ , (?)* 1,+T ') (113) 

Any term on the right hand side of (113) has the generic form (108). All the St < S, so that one can 
invoke the inductive hypothesis (112). The actual degree of divergence of this product is related to the 
naive index of its corresponding constituents as: 

da ( II H(St) ) = \ + d o (constituents of JJ H (s ^ if k out of D sets Si are non-empty . (114) 

i i 

This apparent anomaly stems from the fact that if a set Si is empty, the contribution scales like u° 
rather than u 1 / 2 . This shows that one can discard almost all terms in the sum on the right hand side: in 
fact, only three classes of term survive: 



jj(n 1 ,...,n p ) _ 



M 



EIK#+E II 



i j H,»2 JieJl,j26J2 

ii¥=i2 JiLSJ 2 ={ni,...,n p } (115) 



i j \ r=l ^ ' kjtj 
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The first term survives because all indices rij are attached to a single sub- melon M % . Thus, it generates a 
term containing a factor of _HA ni >-">%>). As a result, it is transferred to the left hand side (where it belongs). 

The second term has naive index: do = —p — | rij. On top of that, all indices are attached to exactly 
two sub-melons Ai n and M 12 . Thus, after summation, it leads to terms with exactly two non-empty sets, 
Si 1 and Si 2 . From eq. (114), the actual degree of divergence of these terms is d a = l + c?o = 1 — V~ \ 52 j n j- 

The third term has naive index: d = —p — | 52j rij + ^ . As all indices are attached to a single sub-melon 
it leads after resummation to a term with actual degree of divergence d a = \ — p — | 52j rij + | , hence it 
contributes also to the leading order divergence. 

All the terms one eliminated from equation (115) fall into one of the following categories: 

f Tl " J ( ft ) f Th-i 1 

- Terms involving at least three different melons hj^^ , and h\Ji 3 . They have naive index at 



least — p — | 52 j n 3 . Hence, their actual degree of divergence is at least | — p — ^52 



3 n i 



Terms involving a 52™=i Cr ) r ^ l M i ^ ^M* ^ an< ^ a "m*i ^ or ^ They have naive index at least 
— p — 1 52 . rij- + | and involve at least two distinct melons. Hence, their actual degree of divergence 



is at least 1 — p — IJ] 



3r „ i i 



2 n 3 1 2 ' 

- Terms involving at least two factors 52"= l ("O^fc^^AI* ^' They have naive index at least — p — 
I n j + 1- Hence, their actual degree of divergence is at least | — p — | 52 j % + 1. 

- Terms involving T". They have either i) naive index at least — p — § 52 j 11 j + f ■ Hence, their actual 
degree of divergence is at least \ ~ V~ \ 52 j n i + fi or **) naive index — p — | 52 j n j + 5 but involve 
at least two distinct melons. Hence, their actual degree of divergence is at least 1 — p— | 52j n j > + §• 

Explicit resummation gives: 

fl-("i.-."-)(l-2D[F0] D - 1 ) = «Dp-l)[ffO] I> - 2 £ H (Si) H (s 2 ) 

Sius 2 =s Sins 2 =0 , , 

(116) 

1 D-i \ - \ - 1 n. 

r 



j;nj>l r—1 



where a term has been transferred to the left hand side as indicated and S\nj is the set S less the element 
rij. Note that the first terms on the right hand side appears only if a partition on S into two nonempty 
sublists Si and S2 is exists, that is \S\ > 2. As shown already, the right hand side scales like nj 
and using equation (99), one has that: 

At this stage, one might like to check that the coefficient of this leading order divergence in (69) does 
not magically vanish. To that end, bounds on this coefficient may be extracted from equation (116) in an 
analogous manner to those derived in [47] and these show that it is strictly greater than zero. □ 
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